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鳝 条 效应 软体 采摘 机 械 手 建 模 与 试验 


25 48 FR , 周 俊 , 张 震 , 江 自 真 , 张 保 华 


(南京 农业 大 学 工学 院 ,210031 南京 ) 


摘 要 :针对 现 有 鳍 条 效应 软体 手指 建 模 分 析 方 法 难以 兼顾 高 精度 与 小 运算 量 的 问题 ,本 研究 以 课 
题 组 设计 的 苹果 采摘 软体 机 械 手 为 对 象 ,在 观察 软体 手指 变形 的 基础 上 ,提出 将 鳍 条 效应 软体 手指 
等 效 为 串 连 匀 链 四 杆 机 构 ; 运 用 虚 功 原理 推导 了 给 定 变形 状态 下 软体 手指 各 处 接触 力 与 驱动 力 珑 
计算 方法 ;提出 用 线 弹 性 捏 段 模 型 描述 变形 恢复 力矩 ,利用 改进 的 粒子 群 算法 求解 各 级 四 杆 机 构 捏 
簧 劲 度 系 数 ;基于 有 限 元 模型 ,结合 二 次 开发 建立 了 给 定 关节 转角 和 充气 压力 下 驱动 力 算 的 BP 神 
经 网 络 静 态 模型 。 搭 建 试验 平台 对 3 个 不 同 大 小 的 仿真 苹果 在 3 个 不 同 高 度 位 置 下 进行 抓 取 力 测 
量 试验 。 试 验 结 果 表 明 ,所 建立 的 力学 模型 计算 抓 取 力 相 对 误差 绝对 值 小 于 8.6% ,与 有 限 元 模型 
精度 相当 ;驱动 力矩 测量 试验 结果 表明 关节 神经 网 络 模型 计算 值 与 测量 值 变 化 趋势 相同 ,相对 误差 
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Modeling and experiment of soft robotic picking 
hand based on fin ray effect 


YI Boyang, ZHOU Jun, ZHANG Zhen,JIANG Zizhen, ZHANG Baohua 


( College of Engineering, Nanjing Agricultural University ,210031 Nanjing, China) 


Abstract :Aiming at the problem that it is difficult for the existing modeling and analysis methods of soft 
finger based on fin ray effect to meet the requirements of high precision and small amounts of calculation, 
we proposed approximately regarding the soft finger as concatenated four-bar mechanisms by observing the 
soft finger deformation. The theoretical model to solve contact force and actuating torque was established 
based on the principle of virtual work when the deformation state is given. The deformation recovery torque 
was modeled using linearly elastic torsional spring , and then an improved particle swarm algorithm was pro- 
grammed to solve the stiffness coefficient of the torsional spring of each phalanx. Based on the Abaqus fi- 
nite element model and combined with the secondary development, the relationship between actuating 
torque and pressure, joint angle was obtained , and then fitted with the aid of BP neural network. Last ,an 


experiment platform was built to measure the grip force of the soft robotic picking hand. The results showed 
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that the absolute relative error of the model was less than 8.6% ,which was equivalent to that of the finite 


element model. The values of actuating torque calculated by BP neural network shared the same trend with 


the values measured and the relative error was under 12.7%. 


Key words: soft robotic hand; fin ray effect; principle of virtual work; finite element analysis; particle 


swarm algorithm 


软体 机 械 手 凭借 良好 的 柔顺 性 和 适应 性 ,在 非 
结构 化 环境 中 可 以 对 形状 复杂 、 大 小 不 一 的 对 象 实 
现 刚性 机 械 手 难以 完成 的 安全 可 靠 抓 取 , 具 有 较 高 
的 研究 价值 "中 。 其 中 以 鱼 鱼 为 灵感 的 鳍 条 效应 软 
体 手指 " 因 其 特殊 的 被 动 包 络 力学 特性 ,上 且 制 作 简 
单 ,扩展 .通用 性 好 ,在 诸多 领域 得 到 了 广泛 的 应 用 。 
HEMMING 等 :5 开发 了 一 种 用 于 温室 环境 采摘 甜 要 
的 机 器 人 ,工作 时 鳍 条 效应 软体 手指 严 持 住 甜 椒 后 ， 
利用 刀片 切断 果 梗 完成 采摘 ;KOLLE 等 "设计 了 一 
种 基于 鳝 条 效应 的 持 物 锯 ,用 于 处 理 血 管 .神经 和 软 
组 织 等 脆弱 的 生物 结构 ;FISCHINGER 等 "将 鳍 条 
效应 软体 机 械 手 作为 看 护 机 器 人 的 末端 执行 器 ,可 
以 完成 对 生活 中 常见 物品 的 拿 取 摆 放 等 动作 ,协助 
行动 不 便 的 用 户 。 

这 些 研 究 通常 仅 将 鳍 条 效应 软体 机 械 手 作 为 其 
整 机 设计 的 一 部 分 ,进行 简单 的 设计 应 用 ,对 机 械 手 
的 抓 取 性 能 并 未 进行 深入 研究 。 在 自动 化 收获 过 程 
中 ,机 械 手 抓 取 力 控制 作为 末端 执行 器 研究 的 关键 
技术 之 一 ,对 实现 无 损 采摘 具有 重要 意义 。SHAN 
等 (中 建立 了 鳍 条 效应 软体 手指 伪 刚 体 模 型 ,可 以 较 
好 地 描述 软体 手指 变形 和 力 , 但 该 模型 依赖 大 量 的 
输入 数据 ,形式 复杂 , 故 主要 用 于 结构 优化 ;针对 常 
见 的 力 传感器 难以 安装 在 软体 手指 上 的 问题 ,EL- 
GENEIDY 等 “和 XU 等 分别 将 传感器 安装 于 机 
械 臂 末端 和 抓 取 对 象 表面 ,分 析 鱼 条 效应 软体 手指 
输出 力 与 驱动 位 移 的 关系 ,但 实际 采摘 中 果实 的 形 
状 尺寸 ,与 机 械 手 的 相对 位 置 都 是 事先 未 知 且 变化 
的 。 某 种 意义 上 鱼 条 效应 软体 手指 和 欠 驱 动手 
指 2 均 为 自由 度数 目 大 于 驱动 数目 的 机 构 , 具 有 类 
似 的 被 动 适应 原理 ,然而 后 者 仍 属于 刚性 机 械 手 范 
畴 ,复杂 度 与 前 者 存在 较 大 差距 。 目 前 有 限 元 方法 
发 展 已 较为 成 熟 , 分 析 参 数 设置 正确 时 其 结果 具有 
较 高 的 可 信 度 于 ?3] ,但 受 较 高 的 计算 成 本 限制 , 难 
以 作为 实时 反馈 控制 的 手段 。 

此 前 的 力学 模型 或 有 较 多 简化 ,忽略 了 抓 取 对 
象 几 何 和 位 置 的 多 样 性 ,或 形式 复杂 计算 量 大 ,不 适 
用 于 实时 控制 。 本 研究 在 此 前 研究 的 基础 上 ,以 课 
题 组 设计 的 鳍 条 效应 软体 采摘 机 械 手 为 对 象 , 观察 
软体 手指 接触 物体 受 力 变形 状况 ,提出 将 软体 手指 


等 效 为 串 连 贸 链 四 杆 机 构 ; 结 合 虚 功 原理 和 改进 的 
粒子 群 算法 建立 了 鳍 条 效应 软体 手指 力学 模型 , 利 
用 有 限 元 模型 和 二 次 开发 ,运用 BP 神经 网 络 对 驱动 
力矩 进行 拟 合 建 模 ;最 后 搭建 试验 平台 测量 采摘 机 
械 手 对 仿真 苹果 的 抓 取 力 以 及 关节 的 驱动 力矩 , 检 
验 所 建立 力学 模型 的 准确 性 。 


1 软体 采摘 机 械 手 设计 


1.1 采摘 机 械 手 整体 结构 


图 1 展示 了 采摘 机 械 手 装配 结构 ,机 械 手 包括 
基 座 气动 弯曲 关节 、 鳍 条 效应 软体 手指 等 部 件 。 工 
作 时 ,气动 讨 曲 关节 和 柔性 波纹 管 充 气 伸 长 ,驱动 上 连 
杆 相 对 下 连 杆 转 过 一 角度 ,带动 软体 手指 绕 弯 曲 关 
节 转 轴 旋 转动 作 , 鱼 条 效应 软体 手指 从 下 向 上 逐渐 
贴 合 至 苹果 表面 形成 被 动 包 络 式 抓 取 , 然 后 通过 采 
摘 机 械 手 旋转 来 扭 断 果 梗 实现 果实 采 梗 的 分 离 ,最 
后 滑 台 气缸 伸 出 ,苹果 从 三 指 形成 的 包 络 空间 底部 
落 和 与 基 座 底部 孔洞 相连 的 收集 管 中 ,完成 收获 。 
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图 1 MUCHOS AC HORE ULT 
Fig.1 Soft robotic picking hand based on FRE 


1.2 软体 手指 结构 与 材料 


基于 鳍 条 效应 原理 ,利用 仿生 学 方法 设计 的 软 
体 手指 采用 “V E” 对 称 结构 ,定义 采摘 时 与 苹果 相 
接触 的 一 面 为 前 面 , 另 一 面 为 后 面 。 软 体 手指 前 面 
后 面 之 间 等 间隔 地 设 有 连接 两 者 的 横梁 。 根 据 苹果 
大 小 ,本 研究 设计 软体 手指 高 90 mm , 底 宽 36 mm ,经 
验 表明 ,软体 手指 横梁 数目 过 少 会 导致 手指 抓 持 能 
力 下 降 以 及 包 络 抓 取 时 接触 面 各 位 置 压 力 分 布 不 
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35] ,过 多 则 会 影响 软体 手指 被 动 包 络 的 性 能 , 故 最 终 
决定 软体 手指 横梁 数 为 9, 其 中 定义 安装 端 即 与 连 
接 块 相 固 连 的 横梁 为 第 1 级 ,未 端 即 最 靠近 指 尖 的 
横梁 为 第 9 级 。 软 体 手 指 采用 硅胶 材料 通过 浇筑 法 
制作 ,按照 国标 GB/T 528 一 92 制作 硅胶 拉 伸 试 样 ， 
使 用 CMT6104 型 电子 万 能 试验 机 进行 单 轴 拉 伸 试 
验 ,选用 二 参数 Mooney-Rivlin Zk 4 i99 2^?" ,对 原 
始 数据 计算 处 理 后 通过 Matlab 曲线 拟 合 工具 箱 拟 合 
得 到 本 构 模型 参数 C, =1.267 MPa, C, =1.329 MPa, 
试验 曲线 和 拟 合 曲线 如 图 2 所 示 。 
20 


一 一 试验 曲线 
18[ — 拟 合 曲线 
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图 2 硅胶 材料 应 力 -应 变 曲线 


Fig.2  Stress-strain curve of silicone material 
2 软体 手指 力学 模型 


2.1 基于 虚 功 原理 的 接触 力 计 算 


用 于 软体 手指 制作 的 硅胶 材料 具有 较 强 的 非 线 
性 , 且 软 体 手 指 自身 结构 复杂 ,结合 了 材料 与 结构 双 
重 非 线性 的 鳍 条 效应 软体 手指 理论 上 具有 无 限 多 的 
自由 度 。 观 察 图 3(a) 软体 手指 与 直径 90 mm 球体 
接触 有 限 元 仿真 Mises 应 力 云图 ,不 难 发 现 软体 手 
指 变形 主要 发 生 在 横梁 两 端 与 前 \ 后面 的 连接 处 ,以 
图 3(b) 所 示 软 体 手指 第 5 级 横梁 处 为 例 ,横梁 两 端 
与 前 .后 面 的 连接 位 置 1 处 应 力 为 5.65 MPa, 位 置 3 
处 为 4.37 MPa, 其 他 位 置 如 2 处 为 1.13 MPa, 位 置 4 
处 为 0.03 MPa。 忽 上 略 横 梁 两 端 与 前 .后 面 的 连接 处 
外 其 他 位 置 的 变形 ,对 软体 手指 进行 伪 刚 性 段 划 分 ， 
得 到 图 4 所 示 的 力学 模型 ,划分 后 n+1 横梁 鳍 条 效 
应 软体 手指 等 效 为 n 级 串 连 贸 链 四 杆 机 构 , 对 这 些 四 
杆 机 构 由 软体 手指 安装 端 到 末端 从 1 到 依次 编号 。 

根据 平面 机 构 自 由 上 度 计算 公式 ,n ZEB eget 
四 杆 机 构 具 有 n+1 个 自由 度 。 如 图 4(a) 所 示 , 设 
第 i 级 四 杆 机 构 前 、 上 、 后 \ 下 杆 长 度 分 别 为 5、 Los 
L; 和 1 , 某 一 时 刻 , 可 以 通过 给 定 第 1 级 下 杆 与 地 面 
Jeff Oa 和 各 级 前 杆 与 地 面 夹 角 bn ,0s,…,0, 共 
n+1 个 变量 来 确定 各 级 后 杆 与 地 面 夹 角 0, ,0,,,…， 


第 和 1 卷 


0, 和 各 级 上 杆 与 地 面 夹 角 ba, 04,…, Oas bario 
(941 为 第 i 级 四 杆 机 构 上 杆 与 地 面 夹 角 , 第 i 级 四 
杆 机 构 上 杆 为 第 i+1 级 四 杆 机 构 下 杆 ,此 处 i=1， 
2,…,n - 1, 如 未 经 特殊 说 明 , 本 研究 默认 i=1， 
2, ,n) 3 

ES 


» 


后 面 
HU / 
横梁 i 


(a) 整体 接触 变形 图 (b) 第 5 级 横梁 处 变形 图 
图 3 Mises 应 力 云 医 


Fig.3 Mises stress contour plot 


(a) 角度 与 杆 长 (b) HRA 
到 4 软体 手指 力学 模型 
Fig.4  Statics model of fin-ray effector 
实际 采摘 过 程 中 , 鳍 条 效应 软体 手指 前 面 与 苹 
果 为 连续 的 接触 ,苹果 对 软体 手指 前 面 的 压强 在 前 
面 与 横梁 连接 的 位 置 处 数值 较 大 ,将 软体 手指 前 面 
连续 分 布 的 压强 离散 化 为 图 4(c) 所 示 的 一 系列 压 
力 太 ,万 , f. ,压力 ff 的 作用 点 为 第 i 级 四 杆 机 构 
前 杆 与 上 杆 的 匀 接 处 ,方向 为 第 i 级 四 杆 机 构 前 杆 
法 线 和 第 i+1 级 四 杆 机 构 前 杆 法 线 的 角 平 分 线 。 
软体 手指 抓 取 苹 果 产 生变 形 后 ,手指 内 部 有 恢复 无 
载荷 时 初始 形态 的 趋势 ,将 软体 手指 等 效 为 串 连 匀 
链 四 杆 机 构 后 ,可 视 为 各 杆 件 受 到 使 机 构 恢 复 初 始 
状态 的 内 力矩 ,这些 内 力 抢 成 对 作用 在 各 级 四 杆 机 
构 前 下 杆 、. 前 上 杆 \ 后 上 杆 和 后 下 杆 之 间 ,两 两 之 间 
大 小 相等 ,方向 相反 ,相互 抵消 。 以 鳍 条 效应 软体 采 
摘 机 械 手 握 拢 为 虚 运 动 正方 向 ,根据 虚 功 原理 ,有 
f'6V 一 T ôb - T 5,66, * Tiu 80, * 
Ti 005, + Tan5bm (1) 


(c) 载荷 与 内 力矩 


其 中 
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fet hh 4. 如 图 5 所 示 , 分析 第 :级 四 杆 机 构 各 杆 虚 角 速 
8V = (8V, 8V, 6V, 8V,)' 度 之 间 的 关系 ,根据 三 心 定 理 , 有 
To =i Ta To Th Te Sbai = (1 + M,)60,, - M50 (3) 
Tru 三 (Ty Tiv Trus Du ôb, = N ôb; + (1 a N;)60,, (4) 

xf el T ds h; 
Ty; (Tiu Teu Teu TA 式 中 . M, ENS ; N, A T sh; 为 第 i 级 四 杆 机 构 
T, = (Tini Ty, Ty, Pu Ty. ) h; 


60m = (805 一 60, ,60p n 60 , :-:,00;, E 60.) 


bru = (805 z 60, ,60,, - 60,,*… ,60,, -60,4) , 
8605, = (80, i5 80,, ,60% z 80,5 ,::,60,, 一 Sba)", 
80,5 m (88, — 80,,,00,, — 60 ,** ,60,, - 60, "uu 


式 中 :了 为 各 级 压力 矢量 f 组 成 的 矩阵 ; 5V 为 各 级 
压力 作用 点 虚 速 度 矢量 OV; 组 成 的 和 矩阵; T, 为 驱动 
JIJE; 60, 为 第 1 级 下 杆 虚 角速度 ; Ti 为 各 级 前 下 
杆 恢复 力矩 Teo: 组 成 的 矩阵 ; Ti 为 各 级 前 上 杆 恢复 
JIJE Tru: 组 成 的 矩阵 ; Too 为 各 级 后 上 杆 恢复 力矩 
Ty; 组 成 的 矩阵 ; Tao 为 各 级 后 下 杆 恢复 力矩 Toy; 组 
成 的 矩阵 ; 604, 为 各 级 前 下 杆 夹 角 变 化 虚 角 速度 
86,5; 组 成 的 矩阵 ; 60, 为 各 级 前 上 杆 夹 角 变 化 虚 角 
速度 Sru: 组 成 的 矩阵 ; 68, 为 各 级 后 上 杆 夹 角 变化 
虚 角 速度 60s; 组 成 的 矩阵 ; Ebo 为 各 级 后 下 杆 夹 角 
变化 虚 角 速度 Sbon: 组 成 的 矩阵 ; 60% 为 第 i 级 前 杆 


虚 角 速度 ; 60,; 为 第 i 级 后 杆 虚 角 速度 ; 60,, 为 第 i 
级 上 杆 虚 角 速度 。 


前 后 杆 的 延长 线 交 点 与 前 上 杆 交 点 的 距离 ; hi 为 第 
i 级 四 杆 机 构 上 下 杆 的 延长 线 交 点 与 前 下 杆 交 点 的 
距离 。 


gr 


图 5 四 杆 机 构 运动 学 模型 


Fig.5 Kinematic model of four-bar mechanism 
根据 正弦 定理 ,有 
Linsin( 60, - 60) 
sin( 60,, — 80; ) 
lasin (80; — 60,,,,) 
= anba — 50, (0) 
将 式 (2) 代 入 式 (1) ,将 式 (5)、(6) 分 别 代 入 式 
(3)、(4) 后 再 代入 式 (1) ,得 
F' = (ToT, $ TuT, + Ty T, 十 TspT,) QE 


h; = 


(5) 


根据 前 文 连续 压强 离散 化 为 压力 的 方法 ,有 压 (7) 
REF =(-T, fi eS 
Jf, — fi tpa qs P Pi] , 则 o VL NM cv 
2 2 80 = (00, Sba 90, Sbp) 
Hess - Xeno, X toostob) (2) T, T, T, T, T; 为 线性 变换 答 阵 ,分 别 为 
-1 1 0 
- (1 +M,) M, 1 0 
T, = - (1 +M,)(1 +M) (1 +M,)M, M, 1 
n-l n-l n-l n-l 
- [1 +m) LII G * 21M, [TI * (21, [TTC + 1n m ol 
i=l iz2 i-3 i-4 
- (1 +M,) 1 4 M, 0 0 E 0 
- (1 +M,)(1 +M) (1 + M)M, 14M, 0 e 0 
T. = - (1 +M,)(1+M,)(1+M,) (1 +M,)(1 +M,)M, (1 + M,)M, 1 + M, zx 0 
- [Ta +m) caen Tarmm [JJO +M)IM, ee 1H, 
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-(N, € Mj) N, +M, 0 0 Wm 0 
- (N, +M,)(1 KM) (N, + M,)M, N, + M, 0 ss 0 
T, = - (N; +M,)(1 € Mj CM) (N; + M) (1 + M,)M, (N, + MM, N, + M, i 0 
- qv, «MO TTG «Mo CN, e MOLTTO «029 (N, + MOLTTG «014 (N, e MOLTTG o 1M, e N, +M, 
- N, N, 0 0 EU 
- N,C1 € Mj) N,M, N, 0 e 0 
T, - NNI +M,)(1 E Mj) N,(1 + M,)M, N,M, N, EU 
n-l n-l n-1 n-l 
- N, [X +M) NUI] CO +M) JM, N, LITT * 21M, NL II +M) IM -- N, 
iel is2 i23 i=4 
1 0 0 0 E 0 
Oa +0 
0 I, cos (727-5 -0,) 0 0 e 0 
7T, =|0 heos( 2 + bs on) I, cos (2275 E A 0 " 0 


经 过 等 效 后 ,根据 式 (7) 即 可 求解 软体 手指 给 
定 变形 状态 下 各 处 的 接触 力 及 需要 的 驱动 力矩 ,为 
了 验证 公式 的 正确 性 ,在 虚拟 样机 软件 ADAMS 中 
建立 如 图 6(a) 所 示 的 串 连 馆 链 四 杆 机 构 模型 ,在 各 
级 四 杆 机 构 杆 件 间 设 置 线 弹性 扭 禾 连接 以 产生 软体 
手指 变形 恢复 力矩 , 扭 管 劲 度 系数 暂 随 意 设置 以 供 
计算 ,本 研究 取 各 捍 策 劲 度 系数 为 0.4 N - mm^, 
动 宁 曲 关节 对 软体 手指 的 驱动 力矩 则 通过 在 第 1 级 
横梁 和 地 面 间 设置 一 负 刚度 扭 答 进 行 模拟 。 在 各 级 
前 杆 施 加 跟随 部 件 旋转 的 压力 ,需要 注意 的 是 ,为 了 
使 压力 大 方向 为 第 ;级 前 村 法 线 和 第 i+1 级 前 杆 法 
线 的 角 平分 线 , 于 ADAMS. 模型 中 将 各 处 压力 拆 角 
为 作用 于 第 ; 级 前 杆 与 上 杆 的 饺 接 处 ,大 小 相等 均 
p DASE ,方向 分 别 为 第 1 级 前 杆 法 线 和 第 i+1 级 
前 杆 法 线 的 两 个 力 , 两 个 力 等 效 的 压力 /; 大 小 为 
Japans cos fra n) , 提交 静态 分 析 运 算 结果 如 
图 6(b) 所 示 。 将 计算 得 到 的 品 连 饮 甸 四 杆 机 构 变 
形状 态 输入 式 (7) 反 求 各 处 接触 力 ,与 ADAMS 软件 
中 对 应 的 等 效 压力 载荷 大 小 一 致 ,驱动 力矩 等 于 负 


刚度 扭 簧 劲 度 系数 绝对 值 乘 以 第 1 级 横梁 转角 ， 
此 可 确定 公式 推导 正确 。 


05, + 05. 05, uu 05, + 05. 
0 I, cos (92 + Srs - on) Ipcos( 2 + frs P on) lacos( 2a + frm : 6.) 


0;, 十 05. 
I cos( 98 + frs Fx 6.) 


(a) ADAMS Tt! 


(b) 计算 结果 
图 6 ADAMS 建 模 验证 
Fig.6 Validation using ADAMS model 


2.2 基于 粒子 群 算 法 的 扭 得 劲 度 系数 优化 


根据 式 (7) , 想 要 计算 软体 手指 工作 时 各 处 的 
接触 力 , 需 要 知道 软体 手指 几何 参数 .变形 状态 和 变 
形 恢复 力矩 。 几 何 参数 和 变形 状态 可 以 通过 观察 测 
量 得 到 ,恢复 力矩 则 与 材料 结构 和 各 位 置 的 变形 大 
小 有 关 ,无 法 直接 测 得 。 根 据 2.1 节 , 各 级 四 杆 机 构 
前 下 杆 .前 上 杆 `\ 后 上 杆 和 后 下 杆 两 两 匀 接 , 且 两 两 
间 存 在 迫使 杆 件 恢复 无 载荷 时 相对 位 置 的 内 力矩 。 
内 力矩 可 以 视 为 各 级 前 下 杆 . 前 上 杆 \ 后 上 杆 和 后 下 
FFEeRERETORCR HIDE TID HG ^ ERS HEIR 
性 可 通过 线 弹性 模型 描述 ,为 了 简化 分 析 , 认 为 同 级 
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四 杆 机 构 中 各 处 扭 纂 劲 度 系 数 相等 。 
ABAQUS 作为 成 熟 的 商用 软件 ,其 建立 的 有 限 
元 模型 具有 较 高 的 可 信 度 ,软体 手指 利用 本 研究 方 
法 和 有 限 元 方法 得 到 的 两 种 模型 ,其 特性 应 该 相近 ， 
本 研究 以 此 为 目标 ,寻求 一 种 扭 纂 劲 度 系数 的 求解 
方法 。 建 立 如 图 7 所 示 的 有 限 元 模型 ,与 2. 1 节 用 
ADAMS 软件 验证 力学 模型 相似 ,该 有 限 元 模型 通过 
施加 力 载荷 得 到 变形 , 设 定 安装 端 第 1 级 横梁 为 完 
全 固定 约束 ,于 软体 手指 前 面 和 每 一 级 横梁 的 连接 
处 施加 垂直 于 前 面 ,跟随 结 点 旋转 的 接触 力 ,提交 运 
算得 到 软体 手指 变形 。7T, T, T T, 均 删 去 第 1 列 
A T, T, T, T, T; 删 去 第 1 行 和 第 1 列 得 到 TT;， 
F WEE 1 行 得 到 了 ,至 此 式 (7) 变 为 
f= (ToT, + ThuTs + T; T, t TjT,) (T) B (8) 
第 ;级 四 杆 机 构 前 下 杆 . 前 上 杆 \. 后 上 杆 、 后 下 
杆 间 的 夹 角 变化 量 分 别 为 Abro 、Abru 、Absu AOp, 
HEIERI K, A 


Tip = (KAO 有 Ag K; Ars), 
Ty; = (RNAb K,AOsw Ks Aga); 
Ti = (KRAO K, AOp K; Abrus)", 
Ty, = (K,A0gy, K, AO, K,A0gs). (9) 


到 7 有 限 元 模型 
Fig.7 ABAQUS model 
给 定 1 组 变形 状态 与 对 应 的 接触 力 ,Ki ,K,,…， 
Ks 为 未 知 量 , 式 (8) 为 包含 8 个 方程 的 线性 方程 组 ， 
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根据 克拉 默 法 则 方程 组 有 且 仅 有 一 组 解 。 然 而 , 仅 
靠 一 组 变形 状态 与 对 应 的 接触 力求 得 的 劲 度 系数 ， 
只 能 使 力学 模型 等 效 该 变形 状态 下 软体 手指 的 力学 
特性 。 为 了 使 力学 模型 能 够 良好 地 模拟 软体 手指 在 
各 种 变形 程度 下 的 力学 特性 ,在 有 限 元 模型 中 施加 
25 组 不 同 的 接触 力 载荷 ,求解 得 到 对 应 的 25 组 变 
形状 态 。8 个 未 知 量 ,200 个 方程 ,对 于 这 样 的 超 定 
方程 组 问题 ,考虑 采用 粒子 群 算法 优化 求解 。 以 式 
(10) , 即 力 学 模型 在 图 8 所 示 25 组 变形 状态 下 计算 
所 得 各 位 置 接触 力 与 对 应 的 有 限 元 模型 接触 力 载 答 
差 值 的 绝对 值 之 和 作为 适应 度 函 数 进行 最 小 值 优 
化 ,这 样 求解 出 的 扭 簧 劲 度 系数 ,可 以 使 力学 模型 与 
有 限 元 模型 更 好 的 相近 。 


110 


100 


90 


80 


70 


60 


Y/mm 


50 


40 


20 Mh, 


X/mm 


8 25 组 变形 状态 
Fig.8 25 deformation states 


Fo = 2 2 f.i) -fiealj,i)| 0) 
RP: fua Q.D 为 第 j 组 变形 状态 有 限 元 模型 施加 
的 第 i 个 接触 力 载荷 ; /(j,i) 为 第 7 组 变形 状态 力学 
模型 求解 的 第 i 个 接触 力 计算 值 。 
粒子 群 算法 具有 简单 易 行 .收敛 速度 快 . 设 置 参 
dubio US ,但 是 在 进行 复杂 高 维 函 数 优化 时 存 
在 着 易 陷入 局 部 最 优 、 迭 代 后 期 粒子 停 沾 和 精度 较 
差 等 缺点 。 针 对 上 述 问题 ,本 研究 在 粒子 群 算法 的 
基础 上 进行 如 下 改进 230 :利用 混沌 序列 产生 大 量 
初始 粒子 ,挑选 优 适 应 度 粒 子 作为 初始 种 群 ; 当 最 优 
粒子 停滞 时 对 最 优 粒 子 进 行 混沌 变换 ;采用 线性 递 
减 惯性 权重 ;将 飞越 边界 的 粒子 随机 置 于 搜索 空间 
内 ,改进 的 粒子 群 算法 流程 如 图 9 所 示 。 
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混沌 序列 产生 初始 粒 村 


y 度 粒子 作为 初始 种 群 


初始 化 粒子 速度 


更 新 粒子 速度 、 位 置 、 计 算 适 应 度 


更 新 粒子 、 种 群 最 优 位 置 


到 9 改进 的 粒子 群 算法 流程 图 
Fig.9 Flow chart of improved PSO algorithm 


2.3 优化 结果 与 分 析 


于 Matlab 环境 编写 算法 程序 ,迭代 得 到 各 级 四 
FOARDIE RZ MAK 10 ,发 现 从 安装 端 到 
末端 ,各 级 四 杆 机 构 扭 簧 劲 度 系数 大 致 呈 减 小 趋势 。 
第 1 级 四 杆 机 构 扭 簧 劲 度 系数 明显 大 于 其 余 各 级 ， 
这 是 因为 软体 手指 安装 端 受 连接 块 约束 ,刚度 大 ;而 
第 4、5 级 劲 度 系数 略微 增 大 ,是 因为 尽管 逐渐 远离 
安装 端 ,但 该 位 置 四 杆 机 构 结构 趋 于 紧凑 ,综合 两 者 
因素 后 刚度 略微 提高 ;而 第 5 级 后 则 靠近 末端 ,与 安 
装 端 距离 过 远 , 所 以 刚度 低 , 劲 度 系数 逐 级 减 小 。 将 
优化 求解 所 得 的 各 级 四 杆 机 构 扭 簧 劲 度 系数 回 代 入 
式 (8) ,根据 式 (11) 可 求 得 各 组 变形 状态 软体 手指 
力学 模型 接触 力 与 有 限 元 模型 各 位 置 接触 力 载 荷 大 
小 之 和 的 相对 误差 yG) 。 如 图 11 所 示 ,本 研究 建 
立 的 力学 模型 计算 接触 力 总 和 相 比 有 限 元 模型 会 偏 
大 ,各 组 相对 误差 均 在 20% 以内。 


S,Mises AVAVAVAVAVAN S,Mises 

(Avg:75%) A (Avg:75%) 
0.00e0 EB 4.70e0 
0.00e0 ES 4.18e0 
0.00e0 ES 3.66e0 
0.00e0 BA 3.15e0 
0.00e0 Ez 2.63e0 

KD 

0.00e0 Eg 2.11e0 
0.00e0 1.59e0 
0.00e0 1.07e0 
0.00e0 7 5.51e-1 
0.00e0 ROON 3 32e-2 


(b) 0.125 MPa 
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8 
2 GG.) - fia G.i)) 
yg) = 一 (11) 
1.0 n 
0.9 
c 08 
- 07 
E 06 
* 05 Es T 
E — 


10 KRIDE BE Zt 


Fig. 10 Stiffness coefficient of springs in four-bar mechanisms 


j 


11 接触 力 相 对 误差 


Fig.11 Relative error of contact force of each group 


2.4 气动 弯曲 关节 静态 特性 建 模 


在 2.1~2.3 节 中 已 建立 鳍 条 效应 软体 手指 的 
力学 模型 ,软体 手指 为 被 动 元 件 ,依靠 气动 弯曲 关节 
驱动 ,为 了 得 到 驱动 力矩 7, 和 充气 压力 了 的 关系 ， 
本 研究 采用 ABAQUS 有 限 元 分 析 软 件 对 关节 进行 
建 模 分 析 , 图 12、 图 13 分 别 展示 了 关节 在 空 载 时 , 即 
驱动 力矩 为 0N . mm 时 ,不 同 充气 压力 下 的 仿真 结 
果 和 关节 在 充气 压力 0.25 MPa 时 ,不 同 转角 bw 下 
的 仿真 结果 。 


S,Mises 

(Avg:75%) 
1.87el 
1.66el 
1.45el 
1.25el 
1.04el 
8.34e0 
6.27e0 
4.21e0 
2.14e0 
7.51e-2 


(c) 0.25 MPa 


图 12 ^ir S T TEASER I F DS CREE FT 


Fig.12 Simulated deformation diagram of pneumatic joint under different pressure 


投稿 网 站 :http://ejam. xjtu. edu. cn. 微 信 公众 号 :应 用 力学 学 报 


S, Mises S, Mises 

(Avg:7596) (Avg:7596) 
1.46e1 1.6lel 
1.29e1 1.43e1 
1.13e1 1.25e1 
9.72e0 1.07el 
8.11e0 8.95e0 
6.50e0 7.17e0 
4.89e0 5.39e0 
3.27e0 3.61e0 
1.66e0 1.83e0 
4.63e-2 5.30e-2 


气动 弯曲 关节 驱动 力矩 大 小 和 充气 压力 关节 
转角 有 关 ,根据 常识 可 知 充 气压 力 一 定时 ,关节 转角 
越 小 驱动 力矩 越 大 。 关 节 转 角 一 定时 ,充气 压力 越 
大 驱动 力矩 越 大 。 为 了 对 三 者 的 关系 进行 量化 分 
析 , 在 关节 有 限 元 模型 的 基础 上 ,利用 ABAQUS 二 
次 开发 接口 Python 编写 程序 ,改变 充气 压力 和 关节 
转角 进行 自动 化 仿真 ,计算 驱动 力矩 ,自动 化 仿真 遍 
Ji 0,0. 01,---,0.25 MPa 共 26 个 气压 和 每 个 气压 下 
的 23 个 转角 ,得 到 598 条 数据 。 以 充气 压力 和 关节 
转角 为 输入 ,驱动 力矩 为 输出 ,搭建 3 层 BP 神经 网 
络 , 将 样本 数据 训练 后 得 到 气动 弯曲 关节 静态 模型 。 
运用 模型 计算 关节 在 给 定 转 角 不 同 气 压 下 的 驱动 力 
^B A RS 14 所 示 , 图 中 共有 5 条 曲线 ,从 下 到 上 
dm H 20° ~ 60" ,以 10° 弟 增 各 转角 下 驱动 力矩 随 
充气 压力 的 变化 关系 ,可 以 发 现 转角 一 定时 驱动 力 
和 矩 随 着 充气 压力 的 增 大 近似 线性 增 大 。 


1 800 


1 500 


1 200 


T,/(N * mm) 
E 


300 
0 
0.025 0.050 0.075 0.100 0.125 0.150 0.175 0.200 0.225 0.250 
P/MPa 
14 不 同 转角 下 驱动 力矩 与 充气 压力 关系 曲线 
Fig.14  Pressure-actuating torque curve of 


joint at different angles 
3 试验 验证 


3.1 机 械 手 抓 取 力 测量 试验 
通过 虚 功 原理 推导 、 粒 子 群 算法 求解 扭 自 劲 度 


(b) 40° 


图 13 气动 弯曲 关节 在 不 同 转角 下 的 仿真 变形 图 


Fig. 13 Simulated deformation diagram of pneumatic joint at different angles 
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S,Mises 

(Avg:75%) 
1.72el 
1.53e1 
1.34e1 
1.15e1 
9.55e0 
7.65e0 
5.7560 
3.85e0 
1.94e0 
4.23e-2 


系数 ,已 经 得 到 了 一 套 完整 的 软体 手指 力学 模型 ,经 
2.3 节 回 代 校 验 ,已 初步 证 明 该 模型 可 以 较 好 地 措 
述 不 同 变形 状态 下 软体 手指 力 与 变形 的 关系 。 为 了 
进一步 验证 模型 的 准确 性 ,制作 软体 采摘 机 械 手 样 
机 并 搭建 试验 平台 (图 15 ) 。 选 取 3 个 形状 较为 标 
WE 大 小 依 序 增 大 的 山东 烟台 红 富 士 苹果 ,进行 三 维 
扫描 ,利用 Geomagic studio 完成 道 向 建 模 设计 ,通过 
3D 打印 得 到 仿真 苹果 ,将 ATI Nano-25E 六 轴 力 / 力 
和 矩 传感器 安装 于 仿真 苹果 内 部 用 于 测量 抓 取 力 。 


图 15 软体 采摘 机 械 手 试验 平台 


Soft robotic picking hand experiment platform 


Fig. 15 

三 指 机 械 手 具有 良好 的 自 定 心 能 力 ,可 以 确保 
苹果 位 于 机 械 手 的 对 称 轴 上 ,但 在 采摘 时 由 于 视觉 
定位 误差 等 原因 ,苹果 相对 机 械 手 的 高 度 位 置 h 会 
在 一 定 范围 内 变化 ,因此 将 3 个 仿真 苹果 依次 通过 
软 绳 悬 吊 于 机 械 手 3 个 气动 弯曲 关节 的 上 、 下 连 杆 
贸 链 所 在 平面 上 方 95 .100 、105 mm ,3 个 不 同 高 度 位 
置 处 ,如 图 16(a) 所 示 , 共 进行 9 次 抓 取 试 验 。 手 动 
调节 调 压 阀 控制 气动 弯曲 关节 充气 压力 ,使 每 次 抓 
取 关 节 转 角 均 为 30°, 图 16(b)、(c)、(d) 展 示 了 机 
械 手 在 3 个 不 同 高 度 位 置 下 对 1 号 苹果 的 抓 取 状 
态 ,可 见 软 体 手指 均 能 对 苹果 实现 良好 的 包 络 贴 合 ， 
25 3 号 苹果 尺寸 大 于 1 号 苹果 ,在 其 余 6 次 试验 
中 软体 手指 也 能 对 苹果 形成 较 好 地 包 络 。 
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例 , 抓 取 力 测量 试验 通过 有 限 元 模型 和 力学 模型 进 
行 仿真 模拟 ,有 限 元 模型 抓 取 力 计算 值 fis 可 以 直 
接 于 后 处 理 模块 读 出 ,力学 模型 抓 取 力 计算 值 FF, 
通过 式 (12) 求 得 。 将 每 组 试验 抓 取 力 测 量 值 ,和 
相对 应 的 Fog, v Fa, 记录 于 图 18 中 。 


Fe 


(12) 

岗 察 试验 结果 ,发 现在 各 条 件 下 鳝 条 效应 软体 
采摘 机 械 手 对 苹果 的 抓 取 力 在 11.33 ~ 19.73 N 之 
间 , 在 同一 高 度 位 置 和 相同 关节 转角 下 ,苹果 大 小 越 
大 ,软体 手指 对 苹果 果实 的 抓 取 力 越 大 ;相同 苹果 在 


(b) 95 mm (c) 100 mm (d) 105 mm 相 同 关 节 转 角 下 , 高 度 位 置 越 大 , 即 越 靠近 软体 手指 
图 16 3 种 高 度 抓 取 力 测量 末端 , 抓 取 力 越 小 。 比 较 各 组 试验 抓 取 力 测量 值 和 
Fig.16 Grip force measurement at three different heights 对 应 的 两 模型 计算 值 ,力学 模型 相对 误差 绝对 值 在 


图 17 以 1 号 苹果 悬挂 在 95 mm 高 度 位 置 处 为 8.6% 以内， ak 


^23 


(a) 实物 试验 (b) 有 限 元 模型 (o 力学 模型 
图 17 抓 取 力 测量 试验 与 仿真 模拟 


Fig.17 Grip force measuring experiment and simulation 
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图 18 ”机械手 对 苹果 在 不 同 高 度 位置 的 抓 取 力 
Fig.18 Grip force at different heights 
采摘 苹果 时 , 机械手 通过 拧 断 实现 果实 与 果 梗 “产生 损伤 , 单 指 抓 取 力 必须 满足 约束 式 
的 分 离 , 拧 断 果 梗 的 扭力 矩 依靠 软体 手指 对 苹果 的 [ < 3 x F x cosa Xu Xr 
摩擦 力 产生 ,考虑 软体 手指 抓 取 力 与 机 械 手 对 称 轴 F < Fan 
的 法 平面 存在 微小 的 夹 角 a 。 要 稳定 采摘 莘 果 上 且 不 RP: 7 为 拧 断 扭力 矩 ; F 为 软体 手指 抓 取 力 ; u 为 


(13) 
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软体 手指 与 苹果 摩擦 系数 ; r 为 苹果 半径 ; Fu NE 
全 抓 取 力 。 


经 摩擦 试验 和 扭转 试验 , 测 得 =0.5, T = 
70 Nmm, 代 入 式 (13) 得 最 小 抓 取 力 为 1. 56 N, 考 
虑 到 一 定 的 可 靠 性 , 抓 取 力 应 不 小 于 5 N; 红 富士 苹 
果 采 摘 期 的 果肉 硬度 P. =1 N/m? ,机械手 抓 取 苹 
果 时 软体 手指 能 与 苹果 果实 表面 充分 接触 ,通过 有 
限 元 后 处 理 模块 读 出 接触 面积 S 约 为 60 mm , 则 最 
大 抓 取 力 即 安全 抓 取 力 f= P xS =60N。 对 比 
试验 结果 ,可 知 软体 采摘 机 械 手 抓 取 力 可 满足 拧 断 
分 离 的 需要 , 且 不 会 对 果实 造成 损伤 。 


3.2 关节 驱动 力矩 测量 试验 


对 气动 弯曲 关节 的 分 析 建 模 是 基于 ABAQUS 
有 限 元 模型 展开 的 ,为 了 验证 分 析 的 可 靠 性 ,进行 关 
节 静 态 试验 ,图 19(a) 为 所 搭建 的 试验 平台 。 调 节 
关节 和 压力 传感器 位 置 ,通过 AS5600 磁 旋 转 编码 器 
模块 保证 关节 转角 为 30°, 从 0. 13 MPa 开始 ,以 
0. 01 MPa 为 步 长 逐渐 增加 充气 压力 至 0.25 MPa, 
次 调节 气压 后 等 竺 气压 传感器 和 压力 传感器 读数 稳 
定 后 记录 关节 驱动 力矩 ,以 避免 柔性 波纹 管 的 迟滞 
性 对 结果 造成 影响 ,试验 结果 如 图 19(b) 所 示 。 


(a) 试验 平台 
1 500 
一 = 一 理论 计算 值 


一 A 一 试验 测量 值 EU 
^ 


N 
[27] 
© 


T,/(N * mm) 
2 
ce 
© 


- 
Cn 
© 


[Iz 
ce 
© 


0.13 0.15 0.17 0.19 021 0.23 0.25 
P/MPa 
(b) 试验 结果 
图 19 气动 弯曲 关节 驱动 力矩 测量 


Fig.19 Actuating torque measurement of joint 
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观察 试验 结果 发 现 ,基于 有 限 元 分 析 的 BP 神经 
网 络 静 态 模 型 计算 的 关节 驱动 力矩 与 充气 压力 呈 较 
好 的 线性 关系 ,模型 计算 值 与 试验 测量 值 整体 趋势 
相同 , 略 大 于 测量 值 , 且 随 着 充气 压力 的 增 大 误差 逐 
渐 增 大 ,原因 在 于 有 限 元 模型 的 材料 .装配 结构 等 无 
法 与 真实 系统 达到 完全 一 致 。 在 0.13 MPa 时 ,相对 
误差 为 最 小 值 4.5% ,在 0. 25 MPa 时 ,相对 误差 达到 
最 大 值 12.7% 。 


4 结论 


1) 基于 二 参数 Mooney-Rivlin 本 构 模 型 对 鳝 条 
效应 软体 手指 进行 仿真 ,发 现 抓 取 时 软体 手指 变形 
主要 发 生 在 横梁 两 端 与 手指 前 后面 的 连接 处 ,提出 
THER T H8 EUN t fe Bet UTILES < 

2) 运 用 虚 功 原理 推导 了 软体 手指 各 处 接触 力 、 
驱动 力矩 与 手指 几何 参数 变形 状态 、 变 形 恢复 力矩 
的 关系 ,以 线 弹 性 扭 簧 模型 描述 恢复 力矩 ,并 用 改进 
的 粒子 群 算法 求解 各 级 四 杆 机 构 扭 簧 劲 度 系数 , 结 
果 表 明 从 安装 端 到 末端 劲 度 系 数 呈 现 减 小 ,略微 增 
大 再 减 小 的 趋势 ;基于 ABAQUS 有 限 元 分 析 、 二 次 
开发 和 BP 神经 网 络 建立 了 气动 弯曲 关节 静态 模型 ， 
得 到 了 关节 驱动 力矩 和 充气 压力 .关节 转角 的 关系 ， 
结果 表明 在 给 定 转角 下 ,驱动 力矩 随 着 充气 压力 增 
大 近似 线性 增 大 。 

3) 搭 建 试验 平台 对 不 同 大 小 苹果 在 不 同 高 度 位 
置 处 进行 抓 取 力 测量 试验 ,并 将 试验 测量 值 与 本 研 
究 所 建立 的 力学 模型 和 ABAQUS 有 限 元 模型 两 种 
模型 的 计算 值 相 比较 。 结 果 表 明 在 同一 高 度 位置 和 
相同 关节 转角 下 ,苹果 越 大 则 抓 取 力 越 大 ,同一 苹果 
在 相同 关节 转角 下 , 莘 果 越 靠近 软体 手指 末端 则 抓 
取 力 越 小 ,软体 手指 抓 取 力 可 满足 无 损 抓 取 并 拧 断 
果 梗 实现 采摘 的 需要 ;本 研究 所 建立 的 力学 模型 计 
算 抓 取 力 相对 误差 绝对 值 不 大 于 8. 6% ,与 有 限 元 
模型 精度 相当 ;测量 30°* 转 角 时 关节 的 驱动 力矩 , 比 
较 发 现 相对 误差 小 于 12.7% 。 

4) 本 研究 建立 的 力学 模型 能 够 根据 变形 计算 鲜 
条 效应 软体 采摘 机 械 手 的 抓 取 力 ,在 保证 精度 的 前 
提 下 , 相 较 于 有 限 元 方法 大 幅 减少 了 运算 量 ,为 鳍 条 
效应 软体 采摘 机 械 手 力 控制 提供 了 理论 基础 。 
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